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We present a case demonstrating the connection between supersymmetric quantum mechanics 
(SUSY-QM), reflectionless scattering, and sohton solutions of integrable partial differential equa- 
tions. We show that the members of a class of reflectionless Hamiltonians, namely, Akulin's Hamil- 
tonians, are connected via supersymmetric chains to a potential-free Hamiltonian, explaining their 
reflectionless nature. While the reflectionless property in question has been mentioned in the liter- 
ature for over two decades, the enabling algebraic mechanism was previously unknown. Our results 
indicate that the multi-solition solutions of the sine-Gordon and nonlinear Schrodinger equations 
can be systematically generated via the supersymmetric chains connecting Akulin's Hamiltonians. 
Our findings also explain a well-known but little-understood effect in laser physics: when a two-level 
atom, initially in the ground state, is subjected to a laser pulse of the form V{t) = (nh/r) / cosh(t/T), 
with n being an integer and r being the pulse duration, it remains in the ground state after the 
pulse has been applied, for any choice of the laser detuning. 

PACS numbers: 46.90.-(-s, 32.80.Qk, 02.30.Jr, 03.65.Fd 



I. INTRODUCTION 

Reflectionless scattering, that is, scattering without 
reflection for all values of the spectral parameter, was 
first studied in the context of hght propagation in a spa- 
tially inhomogeneous dielectric media [1]. Mathemati- 
cally, this is equivalent to finding reflectionless poten- 
tials for the one-dimensional non-relativistic Schrodinger 
equation [3, Isl • The underlying mechanism responsible 
for the reflectionless property of these potentials comes 
from the algebra of supersymmetric quantum mechan- 
ics (SUSY-QM) [J, which links — via a finite number of 
intermediate steps — the reflectionless Hamiltonians to a 
potential- free Hamiltonian . Potential- free Hamilto- 
nians are, in turn, inherently reflectionless. 

It was later discovered that these same reflectionless 
quantum-mechanical potentials, when used as initial con- 
ditions for the Korteweg^ieVries (KdV) equation, lead to 
multi-soliton solutions [l^ . In particular, the reflec- 
tionless potential parameterized by an integer n leads to 
an 71-soliton solution of KdV. 

In this paper we present a SUSY-QM interpretation 
for a known class of reflectionless Hamiltonians [ll|, [13 
that lead to solitons of the sine-Gordon (sG) and attrac- 
tive nonlinear Schrodinger (NLS) equations. We also con- 
jecture a connection between our SUSY chains and a class 
of Darboux transformations [H, [l^l , thus allowing for a 
SUSY-QM explanation for the existence of the latter. 



II. SUPERSYMMETRIC QUANTUM 
MECHANICS 



Imagine we have two Hamiltonians Hq and Hi , which 
are differential operators of finite order and dimension, 
whose coefficients may be functions of space but are 
asymptotically constant as x ^ ±oo, so that a scattering 



problem may be defined for them. Hq and Hi are known 
as supersymmetric partners if they can be factored as 

Ho = AB + e; Hi=BA + e , 

where A and B are also differential operators, and e is 
a constant, multiplied by the appropriate identity. The 
operators A and B also intertwine the Hamiltonians: 

HiB ^ BHq; HqA = AHi . 

Eigenstates of Hq and Hi are related by 

SIV-o) oc \iPi); A\ipi) oc IV-o) , 

where \^o) is an eigenstate oi Hq and j'i/'i) is an eigenstate 
of Hi with the same eigenvalue It is also possible 
to have a hierarchy, or chain, of such supersymmetric 
partners: 
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AiBi+ei 

BiAi + €i = AiBi 
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Hii — BiiAn ~t~ f^n — ■ ■ ■ 

In the case of a SUSY-QM chain, eigenstates of Hn will 
be linked to eigenstates of Hm as 



n<—7n I 'rm I 



(2) 



where T„<_„i = BnBn-i ■ ■ ■ Bm in this particular case. 
More generally, the operator T„2^„j constitutes a so- 
called intertwiner that relates Hamiltonian a Hm to a 
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Hamiltonian _ff„2 . We define an intertwiner as an opera- 
tor that obeys the following property: 

which guarantees that T„2 ^ „j acting on an eigenstate 
of Hn^ produces a linear combination of the eigen- 
states of of the same energy. The intertwin- 

ers can be nested as t„2^„j = t„2^„/T„/^„j 
f ri2 <— n" n" ■(— ri' Tl'n' <— m = In particular: T„^,„ = 
T„^„-iT„_i^„_2 . . . Tm+i<-m- Note that in our case 
T/<-/-i = Bi, hence the expression for T„^„j above. 

Of particular interest is the case in which one of the 
members of the chain, say Hq, is "potential- free" in the 
sense that it is a differential operator whose coefficients 
are constant everywhere. is inherently reflectionless 
since its eigenstates are just plane waves. If the coef- 
ficients in the T;^/_i are asymptotically constant, re- 
lationship ([2|), considered at to = 0, implies that every 
other member of the supersymmetric chain ([1]) is also re- 
fiectionless. In general, Hamiltonians are refiectionless if 
they are connected to a potential-free Hamiltonian via a 
supersymmetric chain. 

III. AKULIN'S HAMILTONIANS 

We studied a class of 2 x 2 matrix differential operators 
given by 

Hn = <yzdx - (TxTi/ cosh(a;) , (4) 
n= -3, -2, -1, 0, +1, +2, +3, ... , 

which we refer to as Akulin's Hamiltonians (see pp. 208- 
210 in [HI). We found a (non-unique) supersymmetric 



Bi+) 

n 

A(+) 

-^71 

with factorization constants 

6(+) = (-l)"(n+i); ,(-) = (-l)"(n-l) . 

Note that different T produce completely different inter- 
mediate Hamiltonians Hn+i/2- The four two- link SUSY- 
QM chains connecting 7J„ to Hn+i via different interme- 



chain linking all the Hamiltonians of form Q. Since Hq 
is potential-free and inherently reflectionless, every other 
member of the chain is also reflectionless. These Hamil- 
tonians are not directly linked, however; there is an in- 
termediate Hamiltonian, which we refer to as i?„+i/2, 

between Hn and Hn+i- The Hamiltonians Hn+i/2 do 
not appear to have any physical significance and further- 
more, they contain spatial-dependencies that diverge as 
X — > ±cx), making the definition of a scattering problem 
impossible. 

We used several symmetries of the Hamiltonians Hn to 
remove some ambiguities in defining such a chain. Con- 
sider the following transformations of operators: Ti = 

O; % = R<Jx- O • (^x^R~\ Ty = (lay) ■ O • 

(icTy)^^; Tz = R(Jz ■ O ■ <^z^R^^^ where R refers to 
the operation x — > —x, and denotes a slot for the 
operator on which the transformation is applied. Un- 
der these actions, the Hamiltonians transform simply: 

Hn ^ +Hn; Hn -Hn] Hn <4 forming 
a one-dimensional representation of the group Dih2 = 
Z2 X ^2, corresponding to 180° rotations about coordi- 
nate axes, plus an identity. As a result, at every step of 
the SUSY-QM chain, the chain quadrifurcates as 

Hn = ... = sr(T[Bi+^]T[M+^]+e^it'>^ 

Hn+1/2 = sr(r[il+)]r[s(+)] + eW) 
Hn+i = (r[Bi;\]r[ii+\] + = ... 

where sj- = T[H]/H, T can be chosen from 
{7}, Tx, Ty, Tz}, and the s values are s/ = -1-1, Sx = -1-1, 
Sy — —1, and Sz = —1. The base SUSY-QM factors are 



I 

diate Hamiltonians Hn+1/2 are illustrated in Fig. [TJ 

The intertwiners linking eigenstates of Hn to eigen- 
states of Hn+i (T„+i4_„), and the intertwiners 
linking eigenstates of Hn to eigenstates of Hn-i 
) can be easily extracted from the chain (O: 



1 — ^n/cosh(x) 

(-l)"cosh(a;) - sinh(x) -dx - ^ ((-l)"(n + 1) 4- ntanh(.x)) 

+dx - I ((-l)"(n-f 1) -f ntanh(a;)) -in/cosh(x) 
(—1)" cosh(a;) — sinh(a::) 



1 



(-l)x 



-dx + \ ((-l)"(n - 1) - ntanh(a:)) +in/cosh(2;) 



(—1)" cosh(a;) + sinh(x) 1 

1 +in/cosh(x) 
(-l)"cosh(a;) -^sinh(a;) <9x + 5 ((-l)"(n- 1) - n tanh(a;)) 
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Ty Hn+i/2Ty 




n+1/2 



FIG. 1: Four SUSY-QM chains for Akulin's Hamiltonians: 
thick arrows correspond to SUSY-QM connections. The 180° 
rotations about the axes OX, OY, and OZ correspond to the 
transformations 71i^ = Rox-Q-Ox^ ^ %/ = {io'x)-Ci-i'io'y)~^ , 
and Tz = Raz ■ O ' o-z^R~^ respectively. 



t„+i^„ = «r (r[Bi+i]'r[ii+']), and t„_i^„ = 

vj- {T[Bl^2^T[An 'j^, where we fix the arbitrary phase 

factor vj- to vj- — T[T]/T, with values vj = +1, 
Vx = —1, Vy — +1, and = —1. Substituting the ex- 
pressions for A's given above, one obtains the fohowing 
expressions for the intertwiners: 



t„+i^„ ^dx-{n+ i)tanh(a;) + ]^{i<Jy)/ cosh(x) 
T„-i^„ =^dx + {n- i)tanh(a;) - i(icry)/ cosh(a;). (6) 



Interestingly, while in the intermediate expressions, the 
intertwiners depend on 7", this dependence disappears in 
the final expressions Thus, like Hamiltonians 77„, 
the intertwiners form a one-dimensional representation 
of the group Dih2. 

This (unexplained) property is nontrivial, because the 
eigenstates of Hn are doubly-degenerate, and thus, gen- 
erally, different versions of the intertwiners can target 
different linear combinations of the eigenstates. 



IV. SOME EXTENSIONS OF THE CLASS OF 
REFLECTIONLESS HAMILTONIANS 



So far, we have shown that Hamiltonians iJ„ = Uzdx — 
cfxWn{x) are reflectionless for Wn{x) = nj cosh(a;), for all 
integer values of n. Here, (Tx,y,z are Pauli matrices. It 
can be easily shown that the class of reflectionless Hamil- 



tonians can be extended to 
1 



c^zdx - -{cr+Wn{x; ^, ?7, xo, 4>) + h.c.) 



where cr+ ^ Ux + io'y- 



nS^ exp[i(?7(a; — xq) + ( 
cosh[^(a; — xq)] 



(7) 



V. INVERSE SCATTERING METHOD 

The inverse scattering method was developed to solve 
the initial value problem for integrable nonlinear partial 
differential equations (NPDE's). Associated with each 
integrable NPDE are two linear differential operators, 
L — L{u) and M — M{u), known as the Lax pair. The 
solution u{x, t) of the NPDE appears as a necessary con- 
dition for L and M to satisfy a Heisenberg-type equation 
L = [M, L] [H, Trivially related to L is a Hamilto- 
nian H that defines a spectral problem = [l6l.ll7|. 
In the case of KdV, L = H . For the sine-Gordon and at- 
tractive nonlinear Schrodinger equations, L — az{H ^ X). 
The eigenvector ip evolves in time through M. 

The procedure for finding u{x,t) is as follows. First, 
find the scattering data a.t t — 0, S{0), for H{u{x,0)). 
Next, use the time evolution through M to find the scat- 
tering data at time t, S{t). Lastly, invert the scatter- 
ing data at time t to find u{x,t). If the H operator 
for a NPDE is reflectionless, then the inversion formula, 
known as the Gel'fand-Levitan-Marchenko integral equa- 
tion, greatly simplifies [H, , reducing to a system of 
algebraic equations, and leads to soliton solutions of the 
NPDE. Thus, finding SUSY-QM chains of refiectionless 
Hamiltonians becomes a method of generating large fami- 
lies of multi-soliton solutions of the corresponding NPDE. 



VI. SOLITONS OF THE SINE-GORDON 
EQUATION 



The sine- Gordon equation is given by 
92 



dxdt 



sin(ii) 



where u = u{x, t), and x and t are light-cone coordinates, 
related to lab-frame coordinates x and thy x — [x + t)/2, 
t = {x — t)/2. Direct-scattering for sG is done via: 



HsG = (^zdx - (^xqix, t ^0) 



(8) 



where q{x,t) = —Ux{x,t)/2. For rj ^ (po = 0, Akulin's 
Hamiltonians (O have this same structure, and they 
give reflectionless direct-scattering initial conditions for 
the sine-Gordon equation. Several known solitonic so- 
lutions of sG can be identified. In particular, Hzpi{£, > 
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u[x',t'] 



-15 -10 



FIG. 2: 2-soliton solution of sine-Gordon equation: scatter- 
ing with Hamiltonian H~2{^ ~ l/'s/S, 0, 0, 0) leads to a two- 
soliton solution where two kinks, initially located at the ori- 
gin, move apart with speeds v — ±1/2 in the stationary lab 
frame, shown above. Other values of ^ yield the same solu- 
tions but viewed from moving frames. 



0, 0, xq, 0) leads to a kink (anti- kink) soliton of initial 
position Xq and velocity v = (1— in the lab 

frame: 

u{x,t) = -4tan"i {exp {t[^{x - xq) + t / (]}) . 

The Hamiltonian -ffip2(^ > 0, 0, xq, 0) leads to a two- 
soliton solution, 



u{x, t) 



4^^^-. ' Tsinh{2[e(x-a;o)+V3e]} \ 
2cosh[-^(a;-xo)+</3C] / 



which consists of two kinks (anti- kinks) that, at < = 
collide at xq with velocities vi = -1-1/2 and V2 ~ —1/2 
in the lab frame (Fig. [2]), and that are observed from a 
reference frame of velocity V — (3^^ — l)/(3^^ -I- 1) with 
respect to the lab frame. 



VII. CONNECTION BETWEEN THE SUSY 
INTERT WINERS AND DARBOUX 
TRANSFORMATIONS 

At this point, it is instructive to compare the inter- 
twiners ([B]) to the (well-known) Darboux transformations 
[T3I [l4| . that also allow the generation of new soliton 
solutions of nonlinear equations from known solutions. 
In the context of the sine-Gordon equation, the Dar- 
boux connection states, in particular, that for any eigen- 
state \ip{X)) of the Hamiltonian Hsg (0) (parametrized 
by some q{x)) corresponding to an eigenvalue A, the state 

\ij{X)) = U\i^{\)) 

is an eigenstate of a new Hamiltonian Hsg (parametrized 
by q{x)) with the same eigenvalue. The intertwining op- 
erator [/ has the form 

^-4a„(C(x) + ^) 



and the new field q is 

g» = -<z(x,<) + Ao(C(x)-l/C(x))) , 

where C,{x) ~ {vi{x) + if2ix))/{<fi{x) - i(p2{x)), and 
{(fi{x),(f2{x))'^ — \(j){Xo)) forms a "fixed" eigenstate of 
Hsg of an eigenvalue Aq. We conjecture that for any 
intertwiner T (|6]) identified above, there exists a Dar- 
boux operator U ([9]), such that T and U coincide on 
the subspace of the eigenstates of Hsg of eigenvalue 
A. For example, the choice of Aq = 1/2 and \4>{Xo)) = 
(exp[-|-a;/2], exp[— 2;/2])'^ leads to the intertwiner 

According to our conjecture, our intertwiners do not 
generate new, previously unknown solutions of the sine- 
Gordon equation. Nevertheless, the intertwiners ([6]) al- 
low one to interpret the Darboux links between various 
Hamiltonians Hsg © as a consequence a hidden SUSY- 
QM structure. Recall that in the case of the KdV equa- 
tion (whose analog of Hsg is the usual stationary lin- 
ear Schrodinger operator), a similar interpretation is well 
known [sj-Q. At the same time, we know of no published 
work which finds a similarly simple explanation for the 
existence of Darboux intertwiners (l9l). 



VIII. SOLITONS OF THE NONLINEAR 
SCHRODINGER EQUATION 

The nonlinear Schrodinger equation with attractive in- 
teractions is given by 

i—u = -^T-^u - 2 u\ u . 
ot ox^ 

Direct-scattering is done with 

HNLS = (^zdx~]j^{cr+q{x,t = {]) + h.c.) , 

with q{x,t) — u(x,t). Akulin's Hamiltonians ([T]), 
H±i{S, > 0, 77, Xq, (j)) produce all onc-soliton solutions: 

exp[ii'] 



u{x, t) = ±^ exp[i(i;a;72 -f 0)] 



coshf: 



where w = 277 is the soliton velocity, is its initial posi- 
tion, (f) is its phase 16;], x' = £^{x — vt — xq), and t = £,^t. 
In turn, H±2{S, > 0, 77, xq, (f>) generates the following 
breather-like two-soliton solutions: 

u{x, t) = ±i exp[i(ra72 + (j))] 

4e**' (cosh(3x') + 36^'*' cosh(x') 



X ■ 



■ 3 cos(8i') -I- 4 cosh(2a;') -I- cosh(4a;') ' 

where v, Xq, and (f) now correspond to the velocity, initial 
position, and the overall phase of the breather, respec- 
tively. All n-soliton solutions of this type were first ob- 
tained, using the inverse scattering method, in Ref. jT^ . 
Note that in the NLS case, as well as in the sG case, mul- 
tisoliton solutions can be generated via Darboux trans- 
formations [U, [l3| . 
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FIG. 3: Sech-shaped laser pulse: our SUSY-QM connection of 
Akulin's Hamiltonians to the potential-free Hq also explains a 
series of specific laser pulse shapes that exhibit a zero ground- 
excited population transfer for any value of the laser detuning. 



IX. THE SECH-SHAPED LASER PULSES 



Akulin's Hamiltonians were first studied by Akulin in 



the context of a two- level, time-dependent system [11 1 
that maps to the spatial scattering problem we have 
considered in this paper: here we present the original 
problem. Consider a two-level atom subjected to a time- 
dependent pulse of the form Veg{t) — V/ cosh (t/r) and 
detuning A. Here V is the amplitude of the pulse, r is 
its duration, and |e) and \g) are the excited and ground 
states, respectively. If we represent the probability am- 
plitudes of the ground and excited states by tjjg and ipe, 
respectively, the dynamics of the system will obey 



'dt' 



^—^9 = 



n 



Tcosh (t/r) 
n , A 
2 



rcosh(i/T) 



It is known that for specific values of the pulse ampli- 
tude, given hy V = nh/r, where n is an integer, the 
transition probability is zero regardless of the detuning 
choice A |TI[ (Fig. [3]). Mathematically, this property was 
not well understood. However, a re-interpretation of this 



problem as finding eigenstates of a reflectionless Akulin's 
Hamiltonian Hn{£. = 1/t, 0, 0, 7r/2) ([7]) with an eigen- 
value A — — iA/2 explains the absence of the population 
transfer in the pulse propagation. 



X. SUMMARY AND OUTLOOK 

We have shown that Akulin's Hamiltonians (HI are 
connected to a potential-free Hamiltonian via super- 
symmetric chains ([5]), explaining their refiectionless na- 
ture. Akulin's Hamiltonians lead to multi-soliton solu- 
tions when used as direct-scattering initial conditions 
for the sine-Gordon and attractive nonlinear Schrodinger 
equations. Additionally, we explain why laser pulses of 
the form V{t) = (nfi/r)/ cosh(t/T) do not transfer pop- 
ulation between the levels of a two-level atom, for any 
choice of detuning. 

The most immediate open question is what specific 
multi-soliton solutions the HnS correspond to for n > 2 
for the sinc-Gordon equation. (The soliton solutions gen- 
erated by Hn for the nonlinear Schrodinger equation are 
given in Ref. Additionally, our analysis only con- 

tains one free-parameter ^, which is not sufficient to fully 
classify the known multi-soliton solutions. Additional 
freedoms may come from different factorization energies 
leading to other SUSY-QM chains (similar to the analy- 
sis in Ref. ^ for the KdV case), and from the inclusion 
of spatial-shifts to the initial conditions at every step of 
the SUSY~QM chains. 
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